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■ Abstract. We present our 3.5PN computation of the (2,2) mode of the gravitational 

QQ ■ wave amplitude emitted by compact binaries, on quasi-circular orbits and in the ab- 

sence of spins. We use the multipolar post-Newtonian wave generation formalism, 
extending and building on previous works which computed the 3PN order This calcu- 
lation required the extension of the multipolar post-Minkowskian algorithm, as well as 
the computation of the mass source quadrupole at 3.5PN order. Our result will allow 
more accurate comparisons to numerical relativity, and is a first step towards the com- 
putation of the full 3.5PN waveform amplitude, which would improve the estimation of 
the source parameters by gravitational wave detectors. 
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■ 1. Introduction 

(N 

O \ Gravitational wave detectors such as VIRGO, LIGO and the planned space-based de- 

tector LISA, need an accurate theoretical prediction for the expected signal of compact 
binary systems. Cross-correlation with a bank of templates will allow on-line detec- 
tion and off-line parameter estimation of the source. This has driven a lot of effort 
for gaining accuracy in our modelling of the signal, both with analytical and numeri- 
^ ■ cal methods. Templates currently in use combine post-Newtonian information for the 

\ inspiral phase (B lanchet 2006) with numerical relativity results for the plunge, merger 

5^ ■ and ringdown phases (|Pretoriusll2005l : lBaker et al.l20d6l : ICampanelli et al.ll2006l), whic h 
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are matched together to provide the full gravitational waveform, e.g. iPan et all (120111) . 

In the current state of the art of the post-Newtonian framework, the equations 
of motion and phase evolution of the binary ar e known up to the 3 .5PN order, while 
the full waveform is known up to the 3PN order (Blan chet et al. P2008) — with as usual 
IPN corresponding to 1 /c^. Extending our knowledge of the waveform up to the 3.5PN 
order would improve the parameter estimation of the source. In p articular , for s uper- 
massive black holes in the LISA range, it has been shown by lArun et all (l2007t ) and 



Trias & SintesI (i2008, ) that a significant improvement on the angular resolution and dis- 



tance measurement was obtained by including higher-order harmonics in the signal. 
Another interesting outcome of higher-order computations will be to compare the re- 
sults obtained with the ones of the num erical relativity co des. 

In this work which is a summary of lFaye et al.l ( 20121) . we address the computation 



of the 3.5PN contribution to the quadrupolar mode (2,2) of the waveform. This mode 
corresponds to the leading order quadrupolar radiation, and it is the mode which is 
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computed the most accurately by numerical codes. Our result constitutes a first step 
towards the computation of the full 3.5PN waveform. 



2. 3.5PN equations of motion in the quasi-circular case 

While addressing the computation of the 3.5PN order of the quadrupolar mode in the 
waveform, we shall need 3.5PN order information on the equations of motion as well, 
which is known from previous works. We model the inspiralling compact binary by 
two point-like particles characterized by their masses mi and mi- Denoting the relative 
position and velocities of the two particles hyx=yi-y2 = rn and v = dxdt, using the 
notation = jc x v/|jir x v| for the unit vector orthogonal to the orbital plane (which is 
fixed in the non-spinning case) and finally completing the triad by /I = f x n, we have 
for quasi-circular orbits 



V 

dv 



rn + rojA, 

• X + ( 



-co 



x + (-+2-)v +O{\0) . 



(1) 

(2) 



Here 0{n) means 0(c'"), and the orbital frequency is oj = (p with ip being the orbital 
phase. We used the facts that r = (9(5), r = (9(10), d> = (9(5) for a quasi-circular orbit 
decaying by the effect of gravitational radiation. 

On one hand, the information about the conservative part of the dynamics for cir- 
cular orbits is contained in the generalized Kepler's law a)(r). This relation has been 
computed in harmo nic coordinates (Blanchet et all l2004l ). in Arnowitt-Deser-Misner 



(ADM) coord inates (iDamour et al.l 20011) . and recently retrieved by effective field the- 
ory methods (IFoffa & Sturanil 120111) . It reads, defining m = m\ + 1112, v = mimi/rn^, 



and the FN parameter y = Gm/rc 
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+ 0(8). (3) 



The constant r'^ is specific to harmonic coordinates and disappears from the physical re- 
sults in the end. On the other hand, the 2.5PN and 3.5PN contributions to the equations 
of motion correspond to radiation reaction which drives the evolution of the orbital 
parameters according to : 
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3. Gravitational waveform for non-spinning binaries 
3.1. Definitions and notations 

Defining an asymptotic Bondi-type coordinate system = (cT, X) dThome [198(f), the 
gravitational waveform may be parametrized by so-called radiative symmetric, trace- 
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free (STF) mass and current multipole moments Ul and Vl according to 
4(7 1 

It \ 

77 7T SahQi ^X)hL-l(jR) \ , (5) 

Ci^v + t J / 

keeping only terms of order 0(1 /R), with R and the radial distance and unit vector, 
and Tr = T - R/c the retarded time. We use a multi-indices notation Al = Ai^,„i,, and 
Pijki is the transverse-traceless projector. Next, the two polarizations of the waveform 
are defined as 



= l(PiPj-QiQj)h^, 



(6) 



hx - \{PiQj + PjQi)hJ^, (7) 

where we take as a convention for the polarization vectors P - and Q = e@, where 
(0, 0) refer to the direction N towards the observer. Decomposing h = h+ -ihx in terms 
of spin-weighted spherical harmonics of weight -2, we define the spherical modes of 
the waveform h^'" by : 

+00 { 

ft- J] 2 /2^'"7!^(0,O). (8) 

£=2 m=-f 

This defines in particular the /i^^ mode of the waveform, the computation of which we 
address in this work. Notice that, from the rotational invariance of the system, it can be 
shown that /i^™ depends on the orbital phase (p as e"™*', so that h^^ corresponds to the 
quadrupolar radiation at twice the orbital frequency. 

3.2. Mode separation for planar binaries 

Now, let us stress an important property, valid for any binary without precession of the 
orbital plane, possibly on an elliptical orbit or including spins aligned or anti-aligned 
with the orbital angular momentum. The parity invariance of the problem and the sym- 
metry of revolution together translate into the following relations for the waveform : 

hJ^(-X,T;-yA,-VA) = hJl(X,T;yA,VA) , (9) 

h(7T-&,<p) = h(&,<p), (10) 

where X is the position of the observer, with angular coordinates (0, 0), and j^, va are 
the positions and velocities of the two massive bodies. We set ^ = ^ - O + n/2, so that 
h depends only on and this relative phase by rotational invariance. 

Using the properties of spin-weighted tensor spherical harmonics, one obtains : 

h^"' oc Ul when £ + m even , 

h^"" oc Vl when £ + m odd. (11) 

In particular, h^^ depends only on Uij and not on Vij, and this dependence reads (with 
1,2, 3 labelling the Cartesian coordinates, and recalling that Uij is STF) : 
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4. Computation of the quadrupole mode 

As we have seen, the computation of the /i^^ mode reduces to the one of the radia- 
tive mass quadrupole . Extending the rnultipolar-post-Minkowskian (MPM) formal- 
ism ("Blan chet & Damoiiij[l986l : lBlanchej|l998l : iBlanchet et aljlloos h to 3.5PN order, 



we express Uij first in terms of canonical moments {Mi,S t), and then of source mo- 
ments {Il,Jl} and gauge moments {Wl,Xl,Yl,Zl) which correspond to higher-order 
corrections. The result of this algorithm reads schematically (the explicit formulas are 
too long to be displayed here) : 

Uij = U'^'' + U]f + Ufj'"^ . (13) 

We find that no cubic interactions between momenta arise at this order, for the radiative 
mass quadrupole. J/J"**' contains the leading order contribution, which is simply the 

second time derivative M.j , as well as higher-order instantaneous interactions between 
the Ml and S l- The tail and memory contributions U]J^ and UfJ'^"^ are hereditary, 
and take the form of integrals over the whole past of the source, with a logarithmic 
kernel in the case of tails. Then, the canonical moments {Ml,Sl} can be reexpressed 
in terms of the source and gauge moments {Il,Jl, Wl,Xl, Yl,Zl}, which themselves 
admit expressions as integrals, over space, of the stress-energy pseudotensor which 
includes both matter and non-linearities in the metric perturbation. At leading order, 
canonical and source moments are identical, and they start differ at 2.5PN order. 

Along with lower-order expressions for the other moments, the quadrupole source 
moment Iij has to be computed to the 3.5PN order. For quasi-circular orbits, it admits 
the following structure, with brackets for STF projection : 

Iij = vm\:R x<'x^> + — S v<' v-'^ + ^ — Crr x^'v^^ + O (8) . (14) 

\ 1 c^r I 

The coefficients Ji and S correspond to the conservative part and are already known 
up to 3PN order (" Blanchet & lyed 12004 ). The 3.5PN contribution, which we had to 



compute, enters only the radiation reaction coefficient Crr, for which we obtain : 

/ 256 1532 \ ^ 

Notice finally that the assumption of a quasi-circular orbit plays a crucial role 
when computing hereditary integrals in the memory and tail terms. 



5. Result 

We absorb logarithmic terms, originating from the tails, by the following change of 
phase variable : 

2GMw, ( aj\ 

ip = ip —In — , (16) 

\a)ol 

which takes in fact the same form as at lower order, but where the ADM mass M 
includes its 2PN corrections; ojq relates to a constant tq introduced in the kernel of the 
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tail integrals. Posing 



IGmvx I6n 



h"" = ^^i^ J—H"" e-'""C , (17) 
/?c2 V 5 



with X = (Gmto/c^)^^^ , we finally obtain the following result for H^^ : 

27 / 107 55 \ ^ 2/ 2173 1069 2047 . 

= l+x\ — — + —v\ + 2nx^'^ + x^\ — — -— — v+ 



42 42 / \ 1512 216 1512 

107;r 347r \ ./27027409 856 4281;: 2jt^ 



c/2/ 107;r 347r \ ./27027409 856 4281;: 
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278185 4l7r2\ 20261 , 114635 , 428 
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33264 96 / 2772 99792 
^,V2,JJI^JJJ^,'_4^\J^«E.'^\A. ,18) 



756 \ 378 162 / \ 27 945 

where je is the Euler's constant appearing at 3PN. The newly computed 3.5PN terms in 
this expression are the ones proportional to x'^^^. We find agr eement for y ^ with the 
test mass limit obtained by black-hole perturbation theory in lTagoshi & Sasaki (1994.) . 

This new result is ready for comparison with high-precision numerical relativ- 
ity outputs. However, the computation of the other modes to the 3.5PN order is not 
straightforward, especially for the (2, 1) mode, and is left for future work. 

Acknowledgments. GF, LB, and BRI thank the Indo-French Collaboration (IFC- 
PAR) under which this work has been carried out. 
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